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?-j ■ Abstract 

Non-Gaussianity may exist in the CDM isocurvature perturbation. We pro- 
vide general expressions for the bispectrum and trispectrum of both adiabatic and 
CN ■ isocurvature pertubations. We apply our result to the QCD axion case, and found 

. a consistency relation between the coefficients of the bispectrum and trispectrum 

, : ~ lO 3 ^ 1 ^] 4 / 3 , if the axion is dominantly produced by quantum fluctua- 

tion. Thus future observations of the trispectrum, as well as the bispectrum, will 
be important for understanding the origin of the CDM and baryon asymmetry. 

o 

On 
O 

> 



1 Introduction 



Recent cosmological observations put limits on deviation from the Gaussian cosmological 
perturbation. The WMAP collaboration reported a constraint on the non-linearity pa- 
rameter as —9 < < HI & t 95% C.L., while improved analysis by another group 
lead to a slightly tight but consistent result, -4 < f^ [) < 80 at 95% C.L. [20 Thus 
at present, the observations are consistent with the Gaussian fluctuation. The standard 
inflation scenario, where the total density perturbation is created by the quantum fluctua- 
tion of the inflaton, predicts negligibly small amount of non-Gaussianity, ~ 0(e,rj), 
where e and r\ are slow-roll parameters which must be much smaller than unity for the 
inflation to last long enough 0, IH El E] • Thus it would be fair to say that the current 
observations are consistent with the standard inflation scenario. 

However, there are scenarios where non-Gaussianity is enhanced and can be observed 
by future experiments. Let us consider a light scalar field, a, which obtains a quantum 
fluctuation during inflation. In the curvaton scenario [7J El [9] , a is responsible for gen- 
eration of the density fluctuation. In this case it is known that non-Gaussianity can be 
significantly large depending on its energy density relative to the total energy density of 
the Universe at the times of its decay PJJ dU [121 EES]- Even in the case where a is not 
responsible for generating the total adiabatic perturbation of the Universe, the observable 
level of non-Gaussianity can be imprinted in the density perturbation |14j . 

Non-Gaussianity may also exist in the isocurvature perturbations [H [151 IB]- This 
possibility was extensively studied recently for the case of the CDM isocurvature [T71 
[TBI [TP"] and baryonic isocurvature perturbations [201 I2T) . This isocurvature type non- 
Gaussianity has distinct effects on the cosmic microwave background (CMB) bispectrum 
[T71 [18] . and the current constraint reads fj^°' = — 5 ± 20 [22] (definition of is 
provided lateid) . For instance, non-Gaussianity of isocurvature type can be generated by 
the axion, a pseudo Nambu-Goldstone boson associated with the spontaneous breakdown 
of Peccei-Quinn (PQ) symmetry, which is introduced in order to solve the strong CP 
problem in the quantum chromodynamics (QCD) [231 [21]- The axion is a light scalar 
field and contributes to the cold dark matter (CDM) of the Universe [22]. In particular, 
the axion can have a large isocurvature perturbation [2H1 [21] • It was found in Ref . [17] 
that large enough non-Gaussianity can be generated without conflicting the bound from 
power spectrum of the isocurvature perturbations. 

In this paper we further investigate non-Gaussianity from isocurvature perturbations, 
including not only the bispectrum but also the trispectrum. So far, the trispectrum of 
the isocurvature perturbation has not been studied. As we will see later, however, the 
trispectrum can become as important as the bispectrum in certain situations. 

This paper is organized as follows. In Sec. [21 a general formalism to calculate bis- 
pectrum and trispectrum including isocurvature perturbations is presented. In Sec. El the 
formalism is applied to the curvaton and axion, and it is shown that both can induce 

1 In this paper we consider only local type non-Gaussianity. 

2 Note that our definition of is same as that used in Refs. [T71 [TB], but differs from that in 
Ref. [22]. 
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large non-Gaussianity. Sec. @]is devoted to discussion and conclusions. 



2 Non-Gausianity from adiabatic and isocurvature 
perturbations 

2.1 Non-linear isocurvature perturbations 

We write the perturbed spacetime metric as 

ds 2 = -M 2 dt 2 + a 2 (t)e^ 7ij {dx l + fidt) (dx j + fidt) , (1) 

where M is the lapse function, $ the shift vector, 7^ the spatial metric, a(t) the back- 
ground scale factor, and if) the curvature perturbation. On sufficiently large spatial scales, 
the curvature perturbation if) on an arbitrary slicing at t = tf is expressed by [28] 

ifj(t f ,x) = N(t f ,U : x) - log^y, (2) 

where the initial slicing at t = is chosen in such a way that the curvature perturbations 
vanish (flat slicing). Here N(tf,ti : x) is the local e-folding number along the worldline 
x = const, from t = ti to t = tf. 

We denote by £ the curvature perturbation if) evaluated on the slice where the total 
energy density is spatially uniform (uniform-density slicing). Let us assume that the 
curvature perturbation is generated by the fluctuation of scalar fields, 5<p a . From Eq. (j2|), 
C can be expanded as 

C ~ NM a + l -N ab 5ct> a 5ct> b + l -N ab M a ^ b ^ c + ^N abcd 5(f) a $<t> b $<f> c $<f> d + ■■■ , (3) 
2 o 24 

where N a = dN/d<fi a and so on, and we have included fourth order terms in 8(f) in order 
to calculate the trispectrum in a consistent way. 

In the similar fashion, we introduce Q to denote the curvature perturbation on a slice 
where pi is uniform (dpi = slicing). Here i = 0,6,7 represents the CDM, baryon and 
photon, respectively. If each component does not exchange its energy with the others, 
Q is known to remain constant for scales larger than the horizon [28J. We define the 
(non-linear) isocurvature perturbation between the z-th and j-th fluid components as 



•S, 3k, g). (4) 

In this paper, we consider the isocurvature perturbation between CDM (baryon) and 
photon, S n (S b7 ). It is useful to define the effective CDM isocurvature perturbation, 

S = S C j + -^-Sjyy, (5) 
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where fi&/f2 c ~ 0.2. It can be expanded as follows, 

S ~ <W + \s ab 8<t> a 8<l> h + l -S ab M a ^ b ^ c + ^- A S aba M a 5<P b ^ c ^ d + ■■■ , (6) 

ZD z4 

up to the fourth order in 5<f). The final curvature perturbation at matter dominated era 
on sufficiently large scales is given by [29J 

C MD = C + \S. (7) 

Here ( is evaluated at deep in the radiation dominated era before the cosmological scales 
enter the horizon, but after the CDM decouples from radiationjf) Notice that effects of 
the adiabatic and isocurvature perturbations on CMB temperature anisotropy are quite 
different. At the Sachs- Wolfe plateau, the CMB temperature fluctuation is given by 
(8T/T) = (1/5)C + (2/5)5*, from which we can see that relative importance of the isocur- 
vature perturbation is enhanced with respect to its effect on the density perturbation. 

For simplicity, we assume that masses of {</> a } are light enough, and that the fluctu- 
ations are Gaussian and independent to each other. Then the correlation functions are 
given by the following form, 

(^4> = (2ir) 3 P s ^k 1 )6(k 1 + k 2 )5 ab (8) 

with 

p s m - f|, o) 

where k denotes a comoving wavenumber, and H in f is the Hubble parameter during infla- 
tion. For later use, we also define the following: 

Aj, s £jy = do) 

Since 5cj) a are Gaussian variables, their three point functions vanish : 

Using the expansion Q and flS]), we can calculate the power spectrum, bispectrum 
and trispectrum of the curvature perturbation (J7j). In the following we derive general 
expressions for the bispectrum and trispectrum including both adiabatic and isocurvature 
perturbations. 

3 The expression (|7|) is valid even for the baryonic isocurvature perturbation. This is because the 
baryon number is conserved and the baryons are non-relativistic. 



3 



2.2 Power spectrum 



First we calculate the power spectrum of the curvature perturbation. We define the 
followings, 

(Q i C^^(2n) 3 5(k 1 + k 2 )P c (k l ), (12) 

(Q i S^)^(2n)- i 6(k 1 + k 2 )P <: s(k 1 ), (13) 

(S k Sj; 2 ) = (27rf8(h + k 2 )P S (h). (14) 

Using and (jSJ), we obtain 

1 r d 3 k' 



I /" drk 

P ( (k) = N a N a P s<t> (k) + -N ab N ab J -—P H {k')P H {\k - k'\), (15) 
P C s(*0 = N a S a P s<j) (k) + -N ab S ab J —P H (k')P H {\k - k'\), (16) 

i r cfik' 

Ps(k) = S a S a P S4> (k) + -S ab S ab J —P H (k')P H {\k - k'\). (17) 



After performing the integration, the spectra P^, P^s and Ps can be expressed as 

P c (k) ~ [iV a iV a + N ab N ab A 2 S(j) HkL)}P S(j> (k), (18) 
P C5 (fc) ~ [iV a 5 a + N ab S ab A 2 Slj) \n(kL)}P Slj} (k), (19) 
P S (A;) ~ [S a S a + S ab S ab A% \n(kL)]P 5(j) (k). (20) 

Here we have introduced an infrared cutoff L, which is taken to be of the order of the 
present Hubble horizon scale [301 El] • The recent WM AP5 result shows that the isocur- 
vature perturbations must be subdominant component and the 95% CL. limits read 
Ps(k*) I P({k*) < 0.19 for the uncorrelated case, and Ps{k*) / P^(k^) < 0.011 for the to- 
tally anti-correlated case, where k* = 0.002 Mpc" 1 is the pivot scale [Tj. Thus we can 
safely take P^(k) ~ N a N a Ps^(k), which is subject to the WMAP normalization condition, 
A 2 ( (Q = (k 3 j2n 2 )P c (K) ~ 2.5 x lO" 9 . 



2.3 Bispectrum 

We define the bispectrum of the curvature perturbation B^ D as 

<cco = + & + **> **)■ (21) 

Notice that this contains four kind of terms taking account of the expression (jTJ), such as 
(CCC); (C(S), ((SS) and (SSS). Let us therefore divide B^ D into four parts depending 
on the origin of each term as 

Bf D (k 1: k 2 , k 3 ) = B ca + B ccs + B css + B sss . (22) 
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Full expression for each term is a bit complicated, and given in the Appendix. We define 
non-linerity parameters as 



%C = lf^[Pdki)Pdh) + 2 perms], 
Bas = tf^lPdkjPdki) + 2 perms], 



5' 

B C ss = h^T 2) [Pdh)Pdk2) + 2 perms], 



(23) 



Bsss = l& 0) [Pdh)Pdk2) + 2 perms], 
and using the results given in the Appendix, we obtain 

lAf ] = (]V ^ )2 [N a N b N ab + (N ab N bc N ca + 2N a N bc N abc )A 2 H ln(k b L)} , (24) 
6 1 

^NL° rl) = , M2 [N a N b S ab + 2N ab N a S b + {3N ab N bc S ca 

+ 2(S a N bc N abc + N a S bc N abc + N a N bc S abc ) }A 2 S(j> ln{k b L) ] , (25) 
6 1 

q/nl = 9^jy jy y l^aS b N ab + 2S ab S a N b + {SS ab S bc N ca 

+ 2(N a S bc S abc + S a N bc S abc + S a S bc N abc ) } A% \n(k b L) ] , (26) 
6 1 

^/nl° = 27(A^ N ) 2 [^ a ^ b ^ ab (S ab S bc S ca + 2S a S bc S abc )A 2 - (f> \n(k b L)] , (27) 

where k b = min{fcj} (i = 1,2,3). The first term in was obtained in Ref. |12j . 

Here we have included terms containing the third derivatives of N and S. In a simple 
case where the scalar potentials are given by the quadratic form, such terms including 
the third derivatives can be neglected. In the case of the standard single field slow-roll 
inflation, Eq. (1241 gives /nl ~ 0(e, rj) where e and rj are slow-roll parameters, as explicitly 
shown in Ref. [221 E], which coincides with earlier investigations [I]. Expressions fl2~4l) - fl2"TI) 
generalize the results of Refs. [TSJ [TS] . 

2.4 Trispectrum 

Similarly to the case of bispectrum, the trispectrum of the curvature perturbation T^ MD 
is defined as 

<CCCO = + h + h + h)T^(k h k 2 , ks, h). (28) 

This contains five kinds of terms, such as (CCCC)> ((((S), (((SS), ((SSS) and (SSSS), 
and hence can be divided as 

Tf D {k 1 , k 2 , k 3 , k 4 ) = + T^ s + T^ss + T^sss + T S sss- (29) 
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Full expressions for these terms are given in the Appendix. Defining non-linerity parame- 
ters tnl and #nl, which characterize the trispectrum including isocurvature perturbations 

as 

T CCCC = r^lP^P^P^h) + 11 perms] 

+ ^t\Pdh)Pdh)Pdh) + 3 perms], (30) 

Teas = 4 c L orl) [P f (A; 1 3)P f (A;3)P f (A; 4 ) + 11 perms] 

+ ^& ll) [Pdh)Pdk2)Pdh) + 3 perms], (31) 
T CC ss = 4 c L or2) [P c (A; 1 3)P c (A;3)Pc(A;4) + 11 perms] 

+ ^? 2) lPdki)Pdh)Pdh) + 3 perms], (32) 

T ( sss = 4 c L ° r3) [P c (A; 1 3)P c (A;3)Pc(^) + H perms] 

54 i ^ (33) 
+ 2^£l r3) [Pdh)Pdk2)Pdh) + 3 perms], 

Tssss = 4l o) [P c (k 13 )P c (k 3 )P ( (h) + 11 perms] 

+ ^Wh)Pdk2)Pdh) + 3 perms], ^ 

where we have used the notation k^ = \ki + kj\. Combining the above expressions with 
the results given in the Appendix, we obtain 



(adi) 1 

t, 



[N a N b N ac N bc + N ab N bc N cd N da A 2 5(f) \n(k t L) 



NL ~ (N a N a f 

+2N a N bd N cd N abc A 2 54> \n(k b L) + (2N a N ab N cd N bcd + N a N b N acd N bcd ) A^ln(fcL)] (35) 

4 C L° rl) = ^N N )3 [ 2 ( N aS b N ac N bc + N a N b N ac S bc ) + AN ab N bc N cd S da A 2 S(/) \n(k t L) 

+ 2(S abc N ad N cd N b + 2N abc S ad N cd N b + N abc N ad N cd S b )A 2 S(f> \n(k b L) 
+ {2(S abc N ab N cd N d + N abc S ab N cd N d + N abc N ab S cd N d + N abc N ab N cd S d ) 
+ (2N abc S abd N c N d + 2N abc N abd N c S d ) } A% \n(kL)] , (36) 

4° r2) = )3 [(N a N b S ac S bc + 2N a N ab S bc S c + 2N a S ab N bc S c + S a S b N ac N bc ) 

+ (4N ab N bc S cd S da + 2N ab S bc N cd S da )A 2 H \n(k t L) 

+ 2{2N abc N ad S cd S b + N abc S ad S cd N b + S abc N ad N cd S b + 2S abc N ab S cd N d )A 2 H \n(k b L) 
+ {2{S abc S ab N cd N d + S abc N ab N cd S d + S abc N ab S cd N d + N abc S ab N cd S d (37) 

+ N abc S ab S cd N d + N abc N ab S cd S d ) 
+ (S abc S abd N c N d + 2N abc S abd N c S d + 2N abc S abd S c N d + N abc N abd S c S d )} A 2 5(f> \n(kL)] , 
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_(cor3) _ 1 



[2(S a N b S ac S bc + S a S b S ac N bc ) + AS ab S bc S cd N da A 2 H hx(k t L) 



NL " 27{N a N a f 

+ 2(N abc S ad S cd S b + 2S abc N ad S cd S b + S a b C S a dS C dNb)Ag ( i ) \D.(kbL) 
+ {2(N abc S ab S C dSd + S abc N ab S c dSd + SafcSahNcdSd + S a b c S a bS cd N d ) 
+ (2S abc N abd S c S d + 2S abc S ab dS c N d )} A% \n(kL)] , (38) 



(iso) 



[S a SbS ac Sbc + SabSbcScdSdaA^ ln(k t L) 



NL " 81(N a N a y 

+2S a Sb d ScdSabcA 2 5<j> ln(fcbL) + i^S a S a bS c dSbcd + S a SbS acd Sb C d) ln(fcL)] , (39) 

for r NL , where fc t = min {kij,ki} = 1,2,3,4), where kf = min{kik 2 } and k 3 = 

mm{kikjki}, and 

+W a N b N cd N abcd A 2 S(j) \n(kL)] , (40) 

S^^^ = 3 (i V 1 iV )3 [ 3N * N b S c N <*° + SabcN a N b N c 

+ 3(S abc N ad N cd Nb + 2N abc S a dN cd N b + iV afec iV ad iV cd ,S 6 ) ln(A; 6 L) 
+3(5 a6cd iV ab iV c iV d + N^SMNa + 2N abcd N ab N c S d )A 2 <t> \n(kL)] , (41) 

!|^T 2) = )3 [3iV a 5 fe 5 c iV afec + 3S afec S a AW c (42) 

+ 3{2N abc N ttd S ed S b + N abc S ad S cd N b + S abc N ad N cd S b + 2S abc N ab S cd N d ) ln(fc 6 L) 
+3(N aicd N ai S c S d + 2N abcd S ab N c S d + 2S abcd N ab N c S d + S abcd S ab N c N d )A 2 H ln(kL)] , 

^fi , NL r3) = 27(iV N ) 3 ^ Sa ^ bN °^ abc + N abcS a S b S c 

+ 3(N a b c S ad S cd Sb + 2S a b c N ad S cd Sb + S a b c S ad S cd Nb) A^ln^L) 

+3y / N abcd S a bS c S d + S a b cd N ab S c S d + 2S a b cd S a bS c N d ) A^ln(fcL)] , (43) 

25#N1? = gl(j\T jy )3 [SaSbScSabc + SSaSbdScdSabcA^ln^kbL) 

+33 a 3 h 3 cd S abcd A%]n{kL)] , (44) 

for (7nl- The first terms in t^l anc ^ #nl were given in Ref. [32], and next order ex- 
pressions were presented in Ref. [H]. In the case of quadratic scalar potentials for <f) a , 
#nl are suppressed compared with tnl since <?nl always include third order derivatives 
such as N abc and S abc . Note that these r NL and #nl are defined through the curvature 
perturbation in the matter dominated era at large scales. As already mentioned, the 
effect on the temperature fluctuation is enhanced at the Sachs- Wolfe plateau in the case 
of the isocurvature perturbation, compared with the adiabatic curvature perturbation. 
For example, at large scales, the primordial non-Gaussianity with = 1 may have 
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about 1000 times enhanced effect on the trispectrum of the temperature fluctuation on 
large scale, than that with = 1. In order to see how they exhibit themselves on the 
CMB trispectrum beyond the Sachs- Wolfe plateau, we must connect the curvature and 
isocurvature perturbations to the temperature fluctuation using the transfer functions, as 
we did in the case of bispectrum [TTJ, [18]. Thus their effects on the CMB temperature 
anisotropy may be quite characteristic. This will be studied in a separate paper. 



3 Application 

We have provided general formulae of /nl, tnl and #nl in the general multi- field case in 
order to keep our arguments most generic. Although these formulae are a bit complicated, 
we see that in this section they can be reduced to a simple form in some models. 

Let us consider a case where both the inflaton <p and another light scalar o obtain 
quantum fluctuations during inflation. We assume that a has a quadratic potential, and 
that the inflaton generated negligible amount of non-Gaussianity. Then terms involving 
Njxj,, N aaa and higher derivatives can be neglected. If the a field decays into radiation 
at some epoch well before the big-bang nucleosynthesis, the primordial fluctuation of a is 
inherited in the adiabatic perturbation, and so is its non-Gaussianity. Then non-linearity 
parameters of the adiabatic type are given by 



6 (adi) NlN aa + iVjA^ ln(k b L) 

5 /nl " m + N*y 



(45) 

and 



(adi ) _ NlNl + N^A% \n(k t L) 

The other non-linearity parameter, ^nl, contains higher derivative terms and hence the 
dominant contribution to the trispectrum comes from tnl- We study this case in the next 
subsection. 

Let us consider another case where a is stable and contributes to the CDM at present 
epoch. Again the potential of a is assumed to be quadratic in a. In this case, the 
fluctuation of a becomes the CDM isocurvature perturbation. Notice that the energy 
density of o is completely negligible when cosmological scales of interest enter the horizon 
in the radiation dominated era, and hence N a , N aa and so on, can be safely neglected. 
Therefore this case results in the non-Gaussianity of isocurvature type, 

6 (iso) _ S 2 a S aa + SlA 2 5<j> ln(k b L) 

g J NL 27 N 4 ' 



and 

C2 C2 I C 
^(iso) _ ^a^aa + *~ 

SUV? 



(iso) _ gg, + S^Aj/, \n(k t L) 
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As we will see below, o can be identified as the axion. In general, calculating the coef- 
ficients S a , S aa , and so on may be a bothersome task. In simple cases, however, we can 
solve the evolution of a analytically (see Refs. [T71 [201 [18] for concrete examples). For 
instance, for the case of quadratic potential, we obtain S a ~ 1/er and S arj ~ 1/cr 2 (cf. 
Eq.([5SD). 

3.1 Cur vat on 

Now let us identify a with a curvaton, which dominantly contributes to the final curvature 
perturbation. We assume that some fraction of the DM (denoted by X) is produced by the 
inflaton decay, and the remaining DM (denoted by Y) originates from thermal bath after 
the curvaton decays. Then X has correlated isocurvature perturbation with the adiabatic 
one, since it fluctuates in the same way as inflaton, not as the curvaton (or equivalently, 
radiation). In this case the fluctuation of a leads to both adiabatic and isocurvature type 
non-Gaussianities. 

For the bispectrum, dominant contributions are /nl^ an d /nl > given as [TTI [T2"| [Tg] 

^Nl di) = ^(3-4/?-2i? 2 ), (49) 

^/nl 0T1) = ~ [3 - AR - 2i? 2 - 2R(1 - e)] (50) 

where R roughly denotes the energy fraction of a field at its decay epoch : R = 3p a / (4p r + 
3p CT ) with p r (p a ) denoting the energy density of the radiation (curvaton) at the epoch 
of curvaton decay, and e denotes the fraction of X to the total DM abundance. From 
WMAP5 constraint on the isocurvature fraction, it is limited as e < 0.035. Here we 
have assumed that classical displacement of the curvaton field is much larger than the 
amplitude of the quantum fluctuation, and hence neglected second term in Eq. (J4"T|) . 
Although is suppressed by small coefficient e, it can have visible effect on the 

temperature fluctuation because of the enhancement of the isocurvature perturbation at 
large scale [T7J Q2] . 

Similarly, we can calculate rj^ an d ^nl^- The result is 

= 4^(3-4i?-2i? 2 ) 2 , (51) 
for the adiabatic type as in the usual case. On the other hand, we also obtain 

4 c ° rl ) = --±(3 - AR - 2R 2 ) [3-AR- 2R 2 - 3R(1 - e)] . (52) 

Other types of t nl 's (t^l , t^ t3 \ rj^) contain more powers of e, and hence have only 
subdominant effect on the trispectrum. In the small R limit, we find another consis- 
tency relation t^ t1 ^ = (48/25) /n^Vnl > w hich, if confirmed, will support the curvaton 
scenario. 



(adi) 
'NL 



6*3 

25 



y.(adi) 
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3.2 Axion 



In this section we apply our results to the axion. The axion, a, is a pseudo Nambu- 
Goldstone boson associated with the spontaneous breaking of the PQ symmetry. The 
breaking scale of the PQ symmetry, F a , is constrained from various experiments as well 
as astrophysical and cosmological considerations. From the observation of SN1987A, 
•^o ~ 10 10 GeV is required [53]. On the other hand, the upper bound is provided by 
the cosmological argument. Since the axion obtains a tiny mass after the QCD phase 
transition, it begins to oscillate coherently since then and it behaves as non-relativistic 
matter. The lifetime of the axion is much longer than the present age of the Universe, 
and hence its coherent oscillations contribute to the present dark matter abundance. The 
abundance of the axion relative to the dark matter abundance, r = Q a /Q c , is estimated 
as [M] 

- 1 - 8X10 " 2 (™)"°' 82 (™) 2 (S)- < 53 > 



where we have defined 



max{F a 9, (54) 



Here 9 denotes the initial misalignment angle of the axion. Thus the PQ scale is con- 
strained as F a < e-^lO 12 GeV. 

If the PQ symmetry is already broken before or during inflation and if it is never 
restored after inflation, the axion obtains quantum fluctuations during inflation and its 
fluctuation becomes the CDM isocurvature fluctuation. This is realized if the inflation 
scale and the subsequent reheating temperature are lower than the PQ scale. Otherwise, 
the axion does not appear during the inflationary era and no CDM isocurvature pertur- 
bation arises. We assume that the inflaton itself does not generate non-Gaussianity, and 
that only the axion has an isocurvature fluctuation. The CDM isocurvature perturbation 
S is given as 



S 



2cii5a ( 5a 
+ — 



al \ a. 



(55) 



where a, = F a 9 denotes the classical deviation from the potential minimum. Here the 
second term in Eq. fl55l) corresponds to non-Gaussian fluctuation. Notice that here we 
have assumed that the axion sits near the potential minimum and the potential is well 
approximated by the quadratic one. While this is a reasonable assumption as will become 
clear, it is still possible that the axion starts oscillation from near the top of the potential. 
In this case the prediction of non-Gaussianity will change [53]. From Eq. (|47p . the non- 
linearity parameter is given by 

^ 1 ^(^[ ffi ' 2 + ^ ln M' (56) 
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The coefficient of the trispectrum, 



is also calculated from Eq. (1481) as 



(iso) 
'NL 



81JV? 



5 ) t< 



(57) 



Figs. [T] and [2] show contours of / N l = 1 an d 100 for F a = 10 10 GeV (black solid lines) 
and r N L 0) = 10 3 , 10 5 and 10 7 (red dashed lines) on H in{ -9 plane, for F a = 10 10 GeV and 
F a = 10 12 GeV. The gray shaded region is excluded from the isocurvature constraint. 
In the blue region the PQ symmetry may be restored during inflation, and so, neither 
isocurvature fluctuation nor non-Gaussianity will arise. When the quantum fluctuation 
overcomes the classical deviation (a 2 <C A 2 a ), the second term in Eq. (156"]) and floTl) 
dominates. In this limit, we find a consistency relation 



(iso) 
'NL 



p (iso) 
J NL 



-i 4/3 



A~ 2/3 C ~ 10 s 



(iso) 
NL 



4/3 



C. 



(58) 



where C = ln(£;t£)/[ln(/cbL)] 4//3 is constant of order unity. This relation is same as what 
is found in Ref. [H] for the "ungaussiton" scenario. It is not surprising the same relation 
holds when the non-Gaussianity comes from isocurvature perturbation if the quantum 
fluctuation dominates, once we look at the similarity of the expression for /nl and tnl 
between the pure adiabatic and isocurvature cases. As stated before, however, it should 
be noted that their effects on CMB anisotropy are completely different. In particular, 
non-Gaussianity of the isocurvature type at large scales is expected to be significantly 
enhanced, and hence these scenarios can be in principle distinguished. 



4 Conclusions 

The present observations show that the density perturbation of the Universe is nearly 
adiabatic, but small admixture of the isocurvature fluctuation is still allowed. Since the 
CDM (baryon) isocurvature perturbation contains information on how it is produced, the 
detection of the isocurvature perturbation and its non-Gaussianity will certainly give us a 
clue to the origin of the density perturbation and underlying physics beyond the standard 
model. In this paper we have derived the formulae for the bispectrum and trispectrum 
of the isocurvature perturbations. Our formulae can be applied to the case where the 
adiabatic and isocurvature perturbations are correlated. As an explicit example, we have 
shown that the QCD axion can have large isocurvature type non-Gaussianity. Therefore, 
the future CMB surveys may be able to detect the trispectrum and the bispectrum of the 
isocurvature perturbation. Such an analysis was performed in Ref. [36] for the adiabatic 
case and it was found that the Planck survey will have a sensitivity t n ^ > 560 at 2a 
level. This cannot be directly translated into the case of isocurvature trispectrum, since 
the effect on the CMB anisotropy is quite different between the adiabatic and isocurvature 
case. We will pursue this issue as future work. 



11 




H inf [GeV] 

Figure 1: Contours of = 1 and 100 for F a = 10 10 GeV (black solid lines) and 

r NL°' ) = 10 3 , 10 5 and 10 7 (red dashed lines). The gray shaded region is excluded from 
iso curvature constraint. In the blue region the PQ symmetry may be restored during 
inflation, and if this is the case, neither isocurvature fluctuation nor non-Gaussianity will 
arise. 



A. Appendix 



In this appendix we collect the form of bispectrum and trispectrum. 
Bispectrum 

The each term in Eq. (122]) is given by 



B ccc (kx, k 2 , fa) = N a N b N ab [Ps^k^Ps^) + 2 perms. 

d 3 k 



+ N ab N bc N ca 
N a N bc N abc 



2tt) 3 
d 3 k 



(27T) 



P54>{k')P54>{\kx + k'\)P 5(k {\k 2 -k'\) 
Ps<p(k')P 5 4,(\h - k'\)Ps4,(k 2 ) + 5 perms 



(59) 
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e 



10 l 



10" 



10" 



10" 



10" 



10" 



10' 



f NL =100 



10 c 




T NL = 10 



10 



12 



Figure 2: Same as Fig. [fl but for F a = 
excluded from the axion overproduction. 



10" 10" 10 
H inf [GeV] 

= 10 12 GeV. The upper blue region 



> 



0.7) is 



B KS (h, h, h) = (N a N b S ab + 2S a N b N ab ) [P&^h)P H (k 2 ) + 2 perms. 

d 3 k' 



+ 3N ab N bc S c 
1 



(27T)" 



+ 7>(S a N bc N abc + N a S bc N abc + N a N bc S abc ) 



x 



d 3 k' 



Ps^Ps^ilh - k'\)P s<j> {k 2 ) + 5 perms 



B(ss(h, k 2 , h) = (S a S b N ab + 2N a S b S ab ) [Ps ( j > (ki)Ps < j > (k 2 ) + 2 perms. 

d 3 k' 



+ 3N ab S bc S c 



(27T) 



+ -(S a S bc N abc + S a N bc S abc + N a S bc S abc ) 



x 



d 3 k' 



P6<f>(k')Ps<f,(\ki - k'\)P S(f> (k 2 ) + 5 perms 



B S ss(h, h, h) = S a S b S ab [P^(fci)P^(Ar 2 ) + 2 perms. 
d 3 k' 



+ S ab S bc S ca 



(27T) 



+ 



S a S bc S, 



abc 



d 3 k' 



Ps^Ps^lh - k'\)P Sl p(k 2 ) + 5 perms 



(60) 



(61) 



(62) 
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Trispectrum 

The each term in Eq. (|29|) is given by 

T^aih, k 2 , k 3 , k 4 ) = N a N b N ac N bc [Ps4>{kx)P S4> {k 2 )Ps^{k 13 ) + H perms. 
+ N a N b N c N abc [Ps^Ps^Ps^ih) + 3 perms.] 

d 3 p 



+ N ab N bc N cd N da 
N abc N abd N c N d 



+ 



+ 



N abcd N ab N c N d 



+ N abc N ab N cd N d 
+ N abc N ad N cd N b 



(2tt) 3 
d 3 p 

(2+J 3 
d 3 p 

d 3 p 



(2tt) 3 
d 3 p 



Ps<j>(p)Ps<f>(\ki - p\)P St p(\k 2 + p\)P s <j>(\p - fci - fc 3 |) + 2 perms. 
P^(p)P,50(|/ci + A; 3 - pl)P^(A; 1 )P ( 5^(A;2) + 11 perms. 
Ps4>(p)Ps^{\k2 + p\)P6<j>(h)Ps^>(h) + 11 perms. 
P5*(p)P54>(\k 2 +p\)Ps4\h + k 2 \)P s ^(k 4 ) + 11 perms. 



(27T) 



3 P60(p)P^(|A;2 +pl)P^(|fei + fc 4 -p|)P^(A;4) + 11 perms. 

(63) 



T« C s(k u k 2 , k 3 , k 4 ) = 2(N a S b N ac N bc + N a N b N ac S bc ) [Ps^Pg^k^Ps^hs) + 11 perms. 
+ (3N a N b S c N abc + S abc N a N b N c ) [Ps^Pg^Ps^h) + 3 perms.] 
+ 4N ab N bc N cd S da 
d 3 p 



x 



+ -{2N abc S abd N c N d + 2N abc N abd N c S d ) 



x 



x 



(2tt) 3 
d 3 p 



P^^P^O&i -pDPs^dfe + p1)P<^(|p - h - k 3 \) + 2 perms. 



Ps<j>{p)Ps<j>{\h + k 3 - p\)Ps4,{ki)Ps < j > (k 2 ) + 11 perms. 



+ -(S ab cdN ab N c N d + N abcd S ab N c N d + 2N abcd N ab N c S d ) 



Ps^{p)Ps^{\k 2 + p\)Ps<t>(k 3 )P S(f) (k 4 ) + 11 perms. 



+ (S abc N ab N cd N d + N abc S ab N cd N d + N abc N ab S cd N d + N abc N ab N cd S d ) 
d 3 p 



x 



(2tt) 



rPs^Ps^lh+pDPs^lh + k 2 \)P 5<t> (k A ) + 11 perms. 



+ (S abc N ad N cd N b + 2N abc S ad N cd N b + N abc N ad N cd S b 
d 3 p 



(2tt) 



TPs<t>{p)Ps<t>{\k2 +p\)Ps<t>(\h + h -p\)P S( p(k 4 ) + 11 perms. 



(64) 
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T a ss(fa, fa, fa, fa) = (N a N b S ac S bc + 2N a N ab S bc S c + 2N a S ab N bc S c + S a S b N ac N bc ) 
x [PstikjPs^k-dPgtfas) + 11 perms.] 
+ (3N a S b S c N abc + 3S abc S a N b N c ) [Ps^Ps^Ps^fa) + 3 perms.] 
+ (4:N ab N bc S cd S da + 2N ab S bc N cd S da ) 
d 3 p 



x 



+ 2 (S abc S abd N c N d + 2N abc S abd N c S d + 2N abc S abd S c N d + N abc N abd S c S d ) 



x 



x 



(2tt) 3 
bdN, 

d 3 p 
(2^)3 

ib 
d 3 p 

W) 1 



Ps4,{p)Ps4>{\fa - p\)Ps<t,(\fa + p\)P S 4,(\p- fa - fa\) + 2 perms. 



P&s{p)P5d>{ | fa + fa ~ p\)Ps4>{fa)P5<j>{fa) + 11 perms 



+ -(N abcd N ab S c S d + 2N abcd S ab N c S d + 2S abcd N ab N c S d + S abcd S ab N c N d ) 



Ps<j,(p)Ps<t>(\fa + p\)Ps<t,(fa)Ps<t>(fa) + 11 perms. 



+ (S abc S ab N cd N d + S abc N ab N cd S d + S abc N ab S cd N d + N abc S ab N cd S d 
+ N abc S ab S cd N d + N abc N ab S cd S d ) 
d 3 p 



x 



(27T) 



f^(p)A^(|fc2 +p|)^(|fci + fc 2 |)^(fc4) + 11 Perms. 



+ {2N abc N ad S cd S b + N abc S ad S cd N b + S abc N ad N cd S b + 2S abc N ab S cd N d ) 
d 3 p 



x 



(27T) 



rP^tjOP^Gfe +p|)P^(|A;i + fc4 -p1)P<^(&4) + H perms. 



(65) 

and T^sss an d Tssss are obtained by interchanging iV and 5 in the expression of P^CCS 
and T^£, respectively. 
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